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FORMATION AND EVOLUTION OF GALAXY DARK MATTER HALOS AND THEIR SUBSTRUCTURE 
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ABSTRACT 

We use the "Via Lactea" simulation to study the co-evolution of a Milky Way-size ACDM halo 
and its subhalo population. While most of the host halo mass is accreted over the first 6 Gyr in a 
series of major mergers, the physical mass distribution [not M v ; r (z)] remains practically constant since 
z = 1. The same is true in a large sample of ACDM galaxy halos. Subhalo mass loss peaks between 
the turnaround and virialization epochs of a given mass shell, and declines afterwards. 97% of the 
z = 1 subhalos have a surviving bound remnant at the present epoch. The retained mass fraction 
is larger for initially lighter subhalos: satellites with maximum circular velocities V m!iX = 10 km/s 
at z = 1 have today about 40% of their mass back then. At the first pericenter passage a larger 
average mass fraction is lost than during each following orbit. Tides remove mass in substructure 
from the outside in, leading to higher concentrations compared to field halos of the same mass. This 
effect, combined with the earlier formation epoch of the inner satellites, results in strongly increasing 
subhalo concentrations towards the Galactic center. We present individual evolutionary tracks and 
present-day properties of the likely hosts of the dwarf satellites around the Milky Way. The formation 
histories of "field halos" that lie today beyond the Via Lactea host are found to strongly depend on 
the density of their environment. This is caused by tidal mass loss that affects many field halos on 
eccentric orbits. 

Subject headings: cosmology: theory - dark matter - galaxies: dwarfs - galaxies: formation - galaxies: 
halos - methods: numerical 



1. INTRODUCTION 

Cosmological N-body simulations with large dy- 
namic range (i.e. with large numbers of parti- 
cles per virialized object and adequately high force 
and time resolution) make it possible to follow 
the highly non-linear formation of cold dark mat- 
ter (CDM) halos and t h eir s u bstruc t ure in great 
detai l (e.g. IGhigna et ail I1998L l2000t lKlvpinjt_alJ 
19991: IMoore et all ll999L~j200ll iFukushige et all 120041: 



in idealized N-body experiments t hat evolve one satel- 
lite in an externa l potential (e.g. [H avashi e t alJ [20031: 
Dekel et al.1 120031 : iKazantzidis et all 12004 iRead et all 



Kravtsov et all 12004: iDiemand et all 120041: IGao et all 
2004 iGill et"atl 12003 : iReed et all 12005ft . We have re- 



cently completed "Via Lactea", the highest resolution 
simulation to date of CDM substructure. The run was 
completed in 320,000 CPU hours on NASA's Project 
Columbia supercomputer, and follows the formation of 
a Milky Way-size halo with 234 million particles, an or- 
der of magnitude more than achieved previously. The 
present-day properties of the galaxy host and its sub- 
structure were presented in Dicmand, Kuhlcn, & Madau 
(2007, hereafter Paper I). In this second paper we use 
data extracted from all 200 snapshots stored during the 
Via Lactea run to study the mass assembly history of the 
main halo and the sub halo p opulation. 

IGhigna et all (|1998ft and iBullock et all (|2001[ ) have 
noted that subhalos are more concentrated than field ha- 
los. We now have the resolution and statistics to quantify 
this effect, and the large number of snapshots allows us 
to understand its origin. We can follow the evolution 
of massive subhalos at similar or higher resolution as 
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2006ft . but within a "live" host halo forming and evolv- 
ing within the cosmological context. The increased res- 
olution allows more accurate estimates of the fraction 
of subhalos that survive, the mass they retain, the ef- 
fect of tidal stripping on their internal structure. How 
subhalo density profile and concentrations evolve dur- 
ing tidal mass loss? Are subhalos really fully disrupted 
once the tida l radius at pericente r is smaller than their 
scale radius ([Havashi et alJ [2003ft ? How strongly docs 
the heating from tidal shocks reduce inner subhalo den- 
sities? Or do tides in the inn er, shallow part of the host 
lead to subhalo compression (|Dekel et al.ll2003ft ? 

One interesting result that becomes apparent in cosmo- 
logical simulations when tracking halos moving within 
a cluster potential is that many (sub)halos that were 
well within the cluster virial radius r v j r at some earlier 
time can be found today beyond r v j r ([Balogh et al.ll2000l : 
IMoore et al.ll2004lGill et al.ll2005ft . This implies that the 
formation histories of "field" galaxy halos are affected 
by their environment, as many of the halos found in the 
outskirts of larger systems today may have shed mass 
during an earlier pericenter passage because of tidal in- 
teractions. Significant correlations between the forma- 
tion times of ga laxy-size halos of a fixed mass, their clus - 
tering strength ([Sheth fc Tormenll2004 IGao et al.ll2005ft 
and the density of their environment ( Harker et al.ll200rf) 
have indeed been found. In this work, we show that these 
correlations are caused by tidal interactions of (sub)halos 
on extended radial orbits with a more massive neighbor. 

This paper is organized as follows. In §[2]we study the 
mass assembly of the Via Lactea host and its substruc- 
ture population, and introduce physical (non-comoving) 
general definitions of (sub)halo properties like size, mass, 
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Fig. 1. — Evolution of radii tm enclosing a fixed mass versus 
cosmic time or scale factor a. The enclosed mass grows in con- 
stant amounts of 0.3 X 10 12 Mq from bottom to top. Shells are 
numbered from one (inner) to ten (outer). Initially all spheres are 
growing in the physical (non comoving) units used here. Shells 
1 to 6 turn around, collapse and stabilize, while the outermost 
shells are still expanding today. Solid circles: points of maximum 
expansion at the turnaround time tt a . Open squares: time after 
turnaround where r?u first contracts within 20% of the final value. 
These mark the approximate epoch of stabilization. The collapse 
factors r M (t ta )/r M (z = 0) for shells 1 to 6 are 3.29, 2.44, 1.98, 
1.70, 1.51 and 1.36, respectively. Thus shells 1 and 2 collapse by 
more than the factor of 2 derived from spherical top-hat, while 
shells 4, 5, and 6 collapse by a smaller factor. 

concentration, and formation time. Section [3] discusses 
the evolution of subhalo concentrations and abundance 
in fixed-mass shells around the host. In §2] we analyze 
individual and ensemble-averaged evolutionary tracks of 
subhalos and discuss how their density profiles evolve 
during tidal shocks at pcricentcr. We present the histo- 
ries and present-day properties of the likely hosts of the 
dwarf satellites around the Milky Way. Average histo- 
ries of subhalos found in certain regions today and their 
survival from z — 1 to z — are also discussed in this 
section. Finally, § [5] summarizes our conclusions. 

2. FORMATION HISTORIES OF GALAXY HALOS 

The Via Lactea sim ulation was performed with the 
PKDGRAV tree-code (|Stadelll200j) and employed mul- 
tiple mass particle grid in itial conditions gen erated with 
the GRAFICS2 package (|Bertschingerll2001f) . The high 
resolution region was sampled with 234 million particles 
particles of mass 2.1 x 10 4 M Q and evolved with a force 
resolution of 90 pc. It was embedded within a periodic 
box of comoving size 90 Mpc, which was sampled at lower 
resolution to account for the large scale tidal forces. We 
adopted the best-fit cosmol ogical parameters f rom the 
WMAP 3-year data release ()Spergel et al.ll2006D : Cl M = 
0.238, fl A = 0.762, H Q = 73 kms" 1 Mpc" 1 , n = 0.951, 
and us = 0.74. The simulation was centered on an iso- 
lated halo that had no major merger after z = 1.7, which 
makes it plausible that this halo would be a suitable host 
for a. Milky Way-like disk galaxy (e.g. iGovernato et alJ 
l2007t ). More details about the Via Lactea run are given 
in Paper I. Movies, images, and data are available at 
http://www.ucolick.org/~dicmand/vl. The host halo 
mass at z = U is -M200 = 1-77 x 10 M0 within a radius 
of T2oo = 389 kpc (we define r2oo as the radius within 
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Fig. 2. — Fraction of material belonging to shell i at epoch a 
that remains in the same shell today. Shells are same as in Fig. 
[T] numbered from one (inner) to ten (outer). Solid circles: time 
of maximum expansion. Open squares: stabilization epoch. Mass 
mixing generally decreases with time and towards the halo center. 

which the enclosed average density is 200 times the mean 
matter density Qm Pent- Note that M200 and T2oo were 
denoted in Paper I as r vir and Mh a io- We revert here to 
the more standard notation for reasons of clarity). 

Before describing the evolution of substructure, we 
have to address the following issues. When does a halo 
become a subhalo? When does the host form and how 
does it grow? Which regions/ volumes should be used 
for a meaningful comparison of subhalo abundances and 
average properties at different cosmic epochs? The com- 
mon procedure is to define at each epoch a "virial" ra- 
dius r v i r , which depends on the cosmic background den- 
sity at the time, and define subhalos as bound clumps 
within this volume. These definitions are not ideal for 
two reasons. First, halos cross this artificial bound- 
ary not only inward ("accretion") but about as often 
also outwards. Averaging over six, relaxed galaxy clus- 
ters, no net infall of subhalos in to the virial region was 
found in iDiemand et al.l (|2004h . and half of the halos 
found today between r v ; r and 2 r v ; r had ac tually passed 
through the cluster at some earlier time ( Baloeh et alJ 
[20001: lMoore"et^[200l iGill et al.ll2005T ). iMaccio et all 
(|2003h noted that the virialized regions of halos are of- 
ten larger than r v ; r , and a lack of mass infall out to 2-3 
virial radii was f ound for a large, rep resentative sample of 
galaxy halos by iPrada et alJ |2006f ). This leads directly 
to the second problem. As the cosmic background den- 
sity decreases with Hubble expansion, formal virial radii 
and masses grow with cosmic time even for stationary 
halos. Studying the transformation of halo properties 
within r v j r (or some fraction of it) mixes real physical 
change with apparent evolutionary effects caused by the 
growing radial window, and makes it hard to disentangle 
between the two. 

To address the second problem, we describe here the 
formation of Via Lactea using radial shells enclosing a 
fixed mass, rjyr- Unlike r v i r , tm stops growing as soon 
as the mass distribution of the host halo becomes sta- 
tionary on the corresponding scale (see Fig. [I]). The 
first problem, however, remains. Mass and substructure 
are constantly exchanged between these shells, as tm is 
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Fig. 3. — Mass accretion history of Via Lactea. Masses within spheres of fixed physical radii centered on the main progenitor are plotted 
against the cosmological expansion factor a. The thick solid lines correspond to spheres with radii given by the labels on the right. The 
thin solid lines correspond to nine spheres of intermediate radii that are 1.3, 1.6, 2.0, 2.5, 3.2, 4.0, 5.0, 6.3 and 7.9 times larger than the 
next smaller labeled radius. Dashed line: A/200- The halo is assembled during a phase of active merging before a ~ 0.37 (z ~ 1.7) and 
remains practically stationary at later times. 

At z = and for a WMAP 3-year cosmology, this yields 
A = 93. Here we introduce the modified formula, 



not a Lagrangian radius enclosing the same material at 
all times, just the same amount of it. The fraction of 
material belonging to a given shell in the past that still 
remains within the same shell today is shown in Figure 
[2j The mixing is larger before stabilization, presumably 
because of shell crossing during collapse. In the station- 
ary phase the shells still exchange mass because many 
particles are on radial orbits. The mixing is smaller near 
the halo center, where most of the mass is in a dynami- 
cally cold, concentrated component that originated from 
the e arliest forming high-er progenitors ([Diemand et al.l 
120051) . 

2.1. Collapse times and collapse factors 

In spherical top-hat collapse, a shell has no kinetic en- 
ergy at turnaround and virializes at half the turnaround 
radius. The final overdensity relative to the critical den- 
sity at the collapse redshift is A = 187r 2 in the Einstein- 
de Sitter model, modified in a flat Universe with a cosmo- 
logical constant to the fitting formula (Bryan & Norman 
1998) 



200 - 82fi A (z)-39fi A (z), 



(3) 



A = 18tt 2 - 82fi A (z) 



39f! A (z), 



where 



n A (z) = 



(i) 



(2) 



and define the virial radius r v j r as the radius enclosing 
a mean density Ap cr ; t . At z = this yields A = 104 
and r v ir = 288 kpc for Via Lactea. We chose this slightly 
different definition for the collapse overdensity so that, 
at high rcdshifts, r v ; r approaches r2oo- 

The simple spherical top-hat collapse ignores shell 
crossing and mixing, triaxiality, angular momentum, ran- 
dom velocities, and large scale tidal forces. Figure Q] 
shows that spheres enclosing a fixed mass have col- 
lapse factors that differ from 2. Inner shells collapse by 
larger factors, in qualitative agreement with the modified 
spherical collapse model of [Sanchez-Conde et al.l (|2006f ) 
that accounts for shell crossing but not angular momen- 
tum. Shells enclosing about the standard viral mass 
collapse by less than a factor of 2, probably because 
of the significant kinetic energy they contain already at 
turnaround. The collapse times are also different from 
spherical top-hat. Shell number five, for example, en- 
closes a mean density of about 104 p cr n today, a virial 
mass of 1.5 x 10 12 M Q and should have virialized just 
now according to spherical top-hat. It did so instead 
much earlier, at a = 0.6. Even the next larger shell 
with 1.8 x 10 12 M Q stabilized b efore a = 0.8 . Our analy- 
sis supports the point made bv lPrada et alj (|2006l ). that 
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spherical top-hat provides only a crude approximation to 
the virialized regions of simulated galaxy halos. 

2.2. Accretion histories 

To understand the mass accretion history of the Via 
Lactea halo we now analyze the evolution of mass within 
fixed physical radii. Figure [3] shows that the mass within 
all radii from the resolution limit of ~1 kpc up to 100 kpc 
grows during a series of major mergers before a = 0.4. 
After this phase of active merging and mass accretion the 
entire system is almost perfectly stationary at all radii. 
The small decrease in density on scales below 1 kpc is 
likely an artifact of time-steps which arc too large com- 
pared to the short dynamical time in these inner regions 
(Paper I). A similar decrease in density in the inner re- 
gions of halos was shown to be caused by too large time 
steps in the convergence tests of Fukushige et al (2004; 
see their Figure 9). Only the outer regions (~ 400 kpc) 
experience a small amount of mass accretion after the 
last major merger, better visible in the linear mass scale 
of Figure |U The mass within 400 kpc increases only 
mildly by a factor of 1.2 from z = 1 to the present. Dur- 
ing the same time the mass within radii of 100 kpc and 
smaller, the peak circular velocity, V max , and the radius 
where it is reached, rymax, ah remain constant to within 
10%. The lack of evolution in the inner density profile, 
and therefore also in V ma , x , r Vm ax and p{< r Vm ax), dur- 
ing this major merger- f ree phase agrees with the findings 
of pr e vious studies fe.g.lWechsler et al.ll2002t I Zhao et all 
1200$ iRomano-Diaz et al.ll2006ft . 

The physical assembly of galaxy halos thus appears 
to occur mainly during an active early phase of ma- 
jor mergers, in which the halo peak circular velocity 
(and the enclosed mass at all radii) grows to its maxi- 
mum value. In contrast to previous wor k on this subject 
(jWechsler et al.ll2002t IZhao et al.|[200l . we do not find 
evidence for much mass growth during the late "slow ac- 
cretion" phase, when the definition of halo mass is based 
on physical instead of comoving scales. Rather, the mass 
distribution and peak circular velocity appear to remain 
constant for the majority of the halo's lifetime. The fact 
that mass definitions inspired by spherical top-hat fail to 
accurately describe the real assembly of galaxy halos is 
clearly seen in Figure |4l where M200 and M v i r are shown 
to increase even when the halo physical mass remains 
the same. This is just an artificial effect caused by the 
growing radial windows r v ; r and r^oo as the background 
density decreases. For Via Lactea M200 increases by a 
factor of 1.8 from z = 1 to the present, while the real 
physical mass within a 400 kpc sphere grows by only a 
factor of 1.2 during the same time interval, and by an 
even smaller factor at smaller radii. 

We find that the small physical accretion since 2 = 1 
seen in Via Lactea is indeed typical of galaxy halos. In a 
45 Mpc periodic box resolved with 300 3 particles of mass 
1.2 x 10 s M (simulated in a 3-year WMAP cosmology), 
we have identified 303 galaxy halos at z — with M200 
ranging from 0.6 x 10 12 M to 5.4 x 10 12 M . The mass 
within a constant physical radius of 400 kpc grows by 
a factor of 1.15+q'j| (the errors indicate the 68% range 
around the median) since z = 1 ("physical accretion"), 
whereas M200 grows by a factor of 2.10+q gg ("apparent 
accretion" ) . For lower mass halos the physical accretion 



is even smaller. In the same 45 Mpc box we find 714 halos 
with M 20 o ranging from 1.5 x 10 11 M to 4.5 x 10 11 M 
today. From z = 1 to the present their mass within 
200 kpc (~ r 2 oo at z = 0) grows by a factor of only 
1.12t°; 2 7, whereas M 20 o increases by 1.85t{5;4o- Physical 
and apparent accretion are correlated, but with a large 
scatter. 

Within the inner 20 kpc, i.e. where the galaxy is ex- 
pected to lie, the gravitational potential remains con- 
stant during the late, quiescent phases of halo formation 
(Figure [3]). Unless there is an evolving, dominant bary- 
onic mass contribution, the rotation rate of the galac- 
tic disk should not evolve in time, with a peak circu- 
lar velocity that may be proportional to the constant 
peak circular velocity of its halo. Assumptions some- 
times made in semi-analytic models about evolution of 
galaxy properties with halo virial quantities (e.g. stellar 
mass with M v - lT ) will produce inaccurate results, consid- 
ering the different length scales and the lack of physical 
accretion both on small and large scales. Models based 
on quantities that do remain constant during stationary 
phases, like peak circular velocity and the corresponding 
radius and enclosed mass, may be more physical. Ob- 
servations of representative samples of z = 1 galaxies, 
including kinematics, are now becomi ng available (e.g. 
IWeiner etlrtl 120061 ; iKassin et "all l2007t ) and it might be 
possible to test whether galaxy radii and masses grow 
like the halo virial scales (i.e. by about a factor of two) 
or if they remain constant like the halo mass distribution 
on physical scales. 

2.3. Formation times 

As discussed above, the common spherical top-hat in- 
spired halo mass definitions M200 and M v - a are not well 
suited to describe the growth of galaxy halos. Care 
should also be used when the complex and extended pro- 
cess of halo formation is quantified with a single number, 
the so called "halo formation time" . Many of the existing 
definitions are based on the evolution of M2 00 or M v - lr 
IWechsler et all 120021: IZhao et all 120031 ; iGao et all 
l2005h . For galaxy halos such formation times depend 
almost exclusively on the amount of apparent accretion 

, which dominates the evolution of M200 and M v ; r for 
more than half of the age of the universe and, in many 
cases, contributes more than half of the total "accreted" 
mass. Since apparent accretion correlates only weakly 
with physical halo growth, it is unclear how and if halo 
formation times calculated in this manner do relate to 
the epoch when most of the physical halo assembly took 
place. Our analysis also casts doubt on whether such for- 
mation times are at all related to the relevant timescalcs 
for galaxy formation. 

Consider, as an example, the widely used definition 
of formation time as the time when M200 (or similarly 
the FOF mass based on a comoving linking length of 0.2 
times the mean particle separations) reaches half of the 
present value. More than half of our large, low-resolution 
sample of galaxy-size halos would form after z = 1 ac- 
cording to this definition, yet their physical mass accre- 
tion is less than 20% over this time span; their mass 

5 The amount of apparent accretion depends on how much mass 
lies in the outer halo, i.e. it is larger for halos with low concentra- 
tions. 
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Fig. 4. — Same as Figure [3] but using a linear scale in enclosed mass. In addition to M200 (upper dashed line) we now also plot M v - lT 
(lower dashed line) and the mass within the radius of maximal circular velocity (dotted line). The physical mass accretion is small after the 
last major merger at a ~ 0.37 (z ~ 1.7): more than 80% of the present-day material within 400 kpc is already in place at z = 1. This value 
is typical for galaxy-size halos. Filled square: median 2 = 1 mass fraction (= 0.87) within 400 kpc for 303 halos of similar mass rescaled to 
today's Via Lactea mass within 400 kpc. Solid circle: corresponding median z = 1 value for M2oo- Filled triangle: corresponding median 
z = 1 value for the mass within 100 kpc. Error bars indicate the 68% scatter around the median. 



assembly was practically completed before their formal 
"formation time". The Via Lactea halo would have a 
formation redshift of z ~ 1 according to this definition 
(see Figure |4} , which is also well after the epoch when 
most of the physical mass accretion actually took place. 
To address this issue, we propose a formation time based 
on peak circular velocity, a quantity that does not evolve 
during the stationary phase of a halo. We define the 
halo formation epoch Zf OTm to be equal to the earliest 
time when V max reaches 85% of its highest value at all 
rcdshifts: 



2.4. A physical (sub)halo concentration index: cy 

Often halo concentrations are presented in terms of 
the virial concentration index defined as the ratio c vu = 
: /r s , where r s is the scale radius of an NFW fit (e.g 



Navarro et al.lll997l; iBullock et all [20011; IWechsler et all 



VmaxOform) = 0.85 max{V max (z)} 



(4) 



Note that a definition of formation time based on the 
present-day peak circular velocity F max (z = 0) would 
lead to significantly higher median formation redshifts, 
since for many halos V max is reduced by tidal stripping. 
We will show in Section POI that this is true even for halos 
beyond the virial radius today, i.e. for "field" halos. For 
comparison, we define the redshift as 



2002HKuhlen et al.H2005t iMaccio' et aljfeoofl ). This def- 
inition has two drawbacks: 1) c vlr grows even during 
epochs of "apparent accretion" , when the physical mass 
distribution remains constant. In Via Lactea, for ex- 
ample, the mass distribution remains nearly unchanged 
from z = 1 to z = 0, but c V i r grows by about a factor 
of two because of the comoving definition of r v ; r ; and 
2) subhalos are truncated at the tidal radius, which is 
always smaller than their formal virial radius, i.e. virial 
radii and thus c v - u - are not well defined for subhalos. 

A direct measure of physical density in the inner re- 
gions of halos is provided b y the "centr a l dens ity param- 



eter" A 



V/2 



is is provided b y the centr a l dens ii 
introduced bv lAlam et~aTl (|2002f) . 



Here we 



refer to this parameter as cy , to avoid confusion with the 
virial overdensity A. 



V max (z 85 ) = 0.85 y max (z = 0) , 



(5) 



but throughout this work we mean Zf or m (eq. HJ when 
we refer to a halo formation time. In Section l4~4l we will 
find a clear environmental dependence in the median z$§ 
of field halos, but not in their median Zform- 



Cy/2 



where r Vn 



p{< fVmax/2) 1 / V m 



Pcrit,0 



2 \-ffo r Vmax/2 



(6) 



nax/2 lS the radius at which the circular ve- 
locity curve reaches half its maximum value. With this 
definition cv/2/°crit,o is equal to the mean physical den- 
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Fig. 5. — Concentration parameters cy (solid) and c v /2 (dotted) 
divided by the density contrast A used to define r v i r at z = 0, as 
a function of c v ; r = r v ; r /r s . 

sity within rvmax/2 and has the advantage that it can be 
directly measured in numerically simulated dark matter 
halos and in observed galactic rotation curves, without 
reference to any particular analytic density profile. Un- 
fortunately, even with Via Lactea's extreme resolution, 
an accurate determination of rvmax/2 is not possible for 
all but the most massive subhalos. Since rvmax is better 
measured, however, we use instead cy, the mean physical 
density within the radius of the peak circular velocity in 
units of /O cr it,0j as our new physical concentration param- 
eter: 



_ P{< IVmax) _ ( V m 
Cy = = Z 



Pcrit,0 



(7) 



For any given analytic density profile it is straight- 
forward to convert b etween c v i r and cy. For an NFW 
([Navarro et al.lll~997l ) density profile, 



p(r) 

the circular velocity is 



r/r s (l + r/r s )' £ 



V c (r)=4irGp s r] 



f(r) 



f(r)=]n(l + r/r s ) 



with 
r/r s 



(8) 

(9) 
(10) 



1 + r/r s 

The maximum of V c (r) occurs at 

n/max = 2.163 r s . (11) 

The NFW scale density p s can be expressed in terms 
of the concentration c v i r and the spherical top-hat virial 
density contrast A (cf. Section l2~Tj) 

^3 



1 ct 



(12) 



-A(z)p crit {z). 
3 /(rvir) 

Combining eqs. ©, (HI), and (O, we find 

c.= (^) 3 4^A(,)^. (13) 
V 2.163/ /(r vir J p crit ,o 

Figure [5] shows a plot of cy divided by the density 
contrast A used to define r v ; r and c v i r , as a function 



of Cvir- For comparison, we also show cy/ 2 /A. Note 
that cy is defined in terms of p cr it today, whereas c v i r 
explicitely depends on redshift through A(z) and p CT it(z). 
The values of cy for a given c v i r that can be read off from 
Figure [5] are thus only valid at z = 0; at higher rcdshifts 
they must be multiplied by [A(z)p cr it(z)]/[A(0)p C rit,o]- 
At z = the Via Lactea host halo has a concentration 
of cy = 3613, which corresponds to c v i r = 10.4. This 
compares well with the value of c v ; r = 11.7 determined 
from the best fitting NFW model. 

3. EVOLUTION OF SUBHALO PROPERTIES 

The present-day cumulative subhalo mass function 
within Via Lactea is well approximated by a simple power 
law 



N(> A/ su b) = 0.0064 



M s 



sub 



-^200 



(14) 



with slope 6 cum — 1 and host halo mass M200 = 1.8 x 
10 12 M Q (Paper I). Here M su b is defined as the mass 

within the tidal radius r t = rcr su b/(v / 2cr ri ost)- This radius 
is the classical Jacobi limit for an isothermal satellite on 
a circula r orbit of radius r within an isothermal host halo 
(see e.g. iRead et al.ll2006t) . It has the property that the 
host local density is half of the local satellite density at 
r t . Similarly the z = subhalo velocity function within 
r2oo is fitted by 



N(> V m ax) = 0.021 



.host 



(15) 



with slope ay ~ 3 and Vmax.host = 182km/s. Here 
we present the time evolution of the normalizations and 
slopes of these two power laws. 

The large number of subhalos in Via Lactea allows us 
to study their abundance, distribution, and concentra- 
tions as a function of distance from the Galactic center. 
For this analysis we use the ten spherical shells plot- 
ted in Figure [fl each at all times containing a mass of 
0.3 x 10 12 M Q and centered on the main galaxy progeni- 
tor. Figure [5] shows the number of subhalos and the sub- 
structure mass fraction of each shell as a function of time. 
The amount of substructure is very closely linked to the 
formation history of its host: in each shell the number of 
subclumps peaks between the epochs of turnaround and 
stabilization. Some time after a shell has stabilized, the 
abundance of subhalos becomes nearly constant. Most 
of the tidal mass loss at a given radius happens early, 
during a relatively short epoch when the corresponding 
shell is near the end of its collapse and approaching sta- 
bilization, and most subhalos are experiencing their first 
pericenter passage (see £|4.5[) . At low redshift, the mass 
in substructure lost from tidal effects or clumps orbit- 
ing out of the shell is approximately replaced by other 
clumps streaming into the shell. The number of satellites 
with Vmax > 5km/s within the final M200 (ie. within 

6 iV(> Mgub) is not a perfect power law: It becomes steeper 
at large masses, due to dynamical friction, and shallower at small 
masses, due to the gradually increasing importance of numerical 
resolution effects. Therefore the best fit slope depends on the mass 
range and the fitting procedure. For M sub > 200 m p we find c*m = 
0.97 ± 0.03. In the same mass range the differential mass function 
dnfJWjjubJ/dJWgub has a<jM = 1-90 ± 0.02. These best fit values 
differ by less than unity, because the differential mass function 
gives more weight to the relatively poorly resolved low mass end. 
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Fig. 6. — Abundance of (sub)halos versus time in shells containing a fixed mass and centered on the main Via Lactea progenitor at 
each time. Shells are ordered from inner (top left) to outer (bottom right). Solid line: number of subhalos with Vmax > 5 km/s. Dashed 
line: number of subhalos with M su ^ > 4.0 X 10 e Mq .Dotted line: mass fraction in resolved (sub)halos within each shell (excluding the 
most massive subhalo to avoid spikes as it orbits though the shells). The vertical dashed line marks the time of maximal expansion of the 
corresponding mass shell, and the vertical solid line the approximate stabilization epoch (see Figure [TJ. The subhalo mass loss rate peaks 
between these two epochs and declines after the region stabilizes. 



all six inner shells) remains within 25 % of its value 
at z = 0.5 at all times until the present. This nearly 
constant substructure abundance over the last 5 Gyrs 
seems to be at odds with the recently found trends that 
old virialized systems are l ess clumpy (iGao et alJl2004t 
Ivan den Bosch et all 120051 ; IZentner et al.l I2005D . How- 
ever, we too see such a trend for the larger subhalos: 
the number of subhalos with K nax > 10 km/s within the 
final M 2 oo decreases steadily from 159 at z — 0.5 to 112 
at z = 0. It seems that the larger subhalos (roughly 
within three decades in mass of the host) are a transient 
population that declines continuously after the host has 
stabilized 7 . On smaller scales, however, substructure 
appears to be more persistent and less dependent of the 
age of the host (cf. iTaffoni et ail (|2003f ) and f4~6]) . 
Today, the number of subhalos with M su b > 4 x 



10 6 Mq is smaller in shells closer to the G alactic center, 
in agreement with p revious studies (e.g. [Ghigna et aTl 
l2lMlDiem"and et all [20041 : IGao et al.ll2004h . The z = 
number density of mass-selected subhalos is well de- 
scribed by a cored isothermal profiles with a scale 
radius similar t o tha t of their host, as proposed by 
iDiemand et al.l (|2004f ) and the ratio of subhalo number 
density and host matter density is simply proportional 
to radius: 

"MoW 



oc r for 0.1 < r/r v j r < 1.0 



(16) 



7 IZentner et al.l 1120051 ) do indeed point out that this trend is 
stronger for more massive subhalos. 



PhostM 

This "spatial antibias" becomes smaller when subhalos 
are selected by their present Vmax, In this case the bias 
scales with enclosed host mass: 

nv ^ ^ oc M(< r) for 0.1 < r/r vil < 1.0 . (17) 
Phost [r) 

The velocity dispersions of these spatially extended sam- 
ples are larger than the dispersions of the more con- 
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Fig. 7. — Evolution of the slopes of the cumulative subhalo velocity (solid lines) and mass (dashed lines) function in the same shells as 
in Fig. [6] Numbers depict the average slopes between the turnaround and stabilization epochs, and from stabilization to the present. The 
slopes show little trend with time or distance from the main progenitor. 

centrated, dark matter component, in good agreement 
with stationary solutions of the Je ans equation applied t o 
this two component system (as in lDiemand et aLll2004j ). 
The orbital anisotropy parameter (3=1 — 0.5of an /of ad 
is identical for both the subhalos and the dark mat- 
ter: (3(r) ~ 0.55(r/r vir ) 1 / 3 for 0.2 < r/r vir < 1.0. 
Tidal mass loss, which causes the spatial and veloc- 
ity distributions of these subhalo samples to differ from 
the dark matter distribution, does not alter the orbital 
properties of substructure. Locally, we find a value of 
P{r = 8kpc) = 0.12. 

While today the distribution of satellites is more ex- 
tended than the dark matter, the opposite is true at 



high redshifts: the inner shells are much more clumpy. 
At a = 0.1 even the smallest halos considered here 
(4.0 x 10 6 Mq) correspond to rare density fluctuations 
(about 2a). The enhanced subhalo abundance in in- 
ner shells at a = 0.1 thus reflects the bias of high-cr 
peaks towards the c e nters of larger scale fl uct uat ions 
(iCole k Kaise71[T989t ISheth fc Tormenlll999h . 

Figure [7] shows the best-fit power-law slopes of the 
substructure cumulative velocity function in the range 
5 — 50 km/s and the cumulative mass function in the 
range 4.0 x 10 6 -4.0 x 10 9 M Q . The innermost shell (top 
left) is affected by numerical resolution effects and small 
subhalo numbers. The slopes show no strong trends with 
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Fig. 8. — Same as Fig. \E\ but now the evolution of the median subhalo concentration (thick line) is plotted versus scale factor. Thin 
line: 68% scatter around the median. All halos with V ma x > 5 km/s are included. 

all radial shells. 8 . In the outer shell the median con- 
centration decreases with time because of the continuous 
formation of new, lower concentration halos. In the inner 
shells the downward trend is halted as the shell collapses 
and the abundance of subhalos freezes (cf. Figure [5]). 
Note how the level of this floor, between turnaround and 
stabilization of a shell, lies at higher concentrations for 
the inner shells as they turnaround at higher redshifts. 
After a shell stabilizes, the median concentration grows 
because tidal forces remove mass from the outer sub- 



time (as in lGao et aLll2004k fReed et al.ll2005| ) or distance 
from the main progenitor. There is a slight trend toward 
steeper mass and velocity functions in the inner shells, 
and toward a steepening with time from the epoch be- 
tween turnaround and stabilization. These small trends 
are suggestive of tidal mass losses being more significant 
in more massive systems and leading to steeper mass and 
velocity functions (see ^4.6)1 . 

The evolution of the median concentration is shown 
in Figure [5J The concentration parameter is related 
to the cosmic mean density at the halo formation time 
(iNavarro et al.lll997t iBullock et aTll2001l ; iWechsler etldl 
l2002f ). Early forming halos, which are well resolved in 
our simulation, do indeed have high concentrations in 



8 With our physical (non-comoving) definition of concentration 
(cq. O, early forming halos have high concentrations from the 
collapse redshift, whereas their c v i r parameter would grow with 
time to reach high values only recently. 
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Fig. 9. — Median subhalo concentrations and 68% scatter (solid 
lines) versus distance from the Galactic center at the present epoch 
(average values over the last ten snapshots from z = 0.05 to z = 0). 
The dotted line shows 400 times the cosmic background density at 
the median formation times of these halos (see Figure fT6t . Finite 
numerical resolution limits concentrations to below a few times 
10 5 p C rit especially in the smaller satellites (the constant upper per- 
centiles within 100 kpc are artificial). Dashed lines are fits (eq. 1186 
to the percentiles measured beyond 100 kpc. Likely dwarf galaxy 
host halos (triangles, same as in Figure [I3t also follow the general 
relation. 

halo regions, thereby reducing rymax and increasing the 
mean density within this radius, cy (jKa zantzidi s et ahl 
12004 see also Section Together with the median, 
the 68% scatter in subhalo concentration is also grow- 
ing. This may be caused by the increasing amount of 
mixing between newly infalling, low-concentration halos 
and strongly stripped, high-concentration clumps. 

At z = we find a clear trend for higher concentra- 
tions closer to the halo center, as shown in Figure [HJ We 
use the following simple empirical fit to approximate this 
relation: 



cv(r) 



Pbgjr) 

Pcrit,0 



(18) 



where pb g is the average density of the corresponding 
spherical shell around the main host. The best-fit coef- 
ficients are (a, b) = (5895, 0.33) for the median concen- 
trations, (2997, 0.16) for the lower (16th) percentile, and 
(19370, 0.39) for the upper (84th) percentile. 

The higher formation rcdshift of the inner halos is not 
enough to fully explain this concentration-radius rela- 
tion. If we simply assume that the median halo concen- 
trations are proportional to the mean cosmic density at 
the median formation epoch of these halos (dotted line in 
Figure 03) we do indeed get a qualitatively correct trend 
with radius, but the effect is not strong enough. Tidal 
interactions must significantly contribute to the final con- 
centration versus radius relation. In ? !4.3l we confirm that 
this is indeed the case: many halos did pass through the 
inner regions of the host at some earlier time and lost 
significant mass from tidal stripping. Interestingly, tides 
seem to increase the median concentrations (and scatter) 



even beyond r2oo = 389 kpc (shells number 7 and 8). To 
summarize, the concentration-radius relation is caused 
by the combined effect of two different processes: 

i) The formation of new small-scale structure stops 
when a shell collapses. Inner shells turn around 
and collapse earlier, and therefore contain earlier 
forming subhalos with a higher median concentra- 
tion. 

ii) Tidal interactions within the host halo increase 
subhalo concentrations. 

4. EVOLUTIONARY TRACKS OF SUBHALOS 



Th e para llel group finder 6DFOF (jDiemand et alJ 
120061 I2007D finds peaks in phase-space density, i.e. it 
links the most bound particles inside the cores of halos 
and subhalos together. The same objects identified by 
6DFOF at different times therefore always have quite a 
large fraction of particles in common. In most cases this 
fraction is over 90% between two subsequent Via Lactea 
snapshots (separated by 68.5 million years). This makes 
finding progenitors or descendants rather easy. When 
tracing halos backwards in time, we link a halo "A" to 
its main progenitor "B" only if A contains at least 50% 
of the particles in B and if B contains at least 50% of 
the particles in A. This definition is time symmetric and 
we use the same links when we follow halo histories for- 
ward in time. When a (sub)halo merges with a larger 
group its forward history ends with a special merger 
flag that points to the ongoing track of the merger rem- 
nant. We include in our analysis only halos larger than 
Knax = 5 km/s at some point during their history. These 
halos are resolved sufficiently well so that one can fol- 
low both their smaller, high-rcdshift progenitors and also 
their present-day remnants, even if they did suffer large 
tidal mass loss. The well resolved sample selected this 
way contains 3883 halos, i.e. it is large enough to offer 
good statistics. Starting at z = we identify the main 
progenitors of all such halos in each snapshot back to 
at least z = 10, when some progenitors start to become 
to small to be resolved and identified with 6DFOF. The 
dotted lines in the right hand panels of Figures [141 and [TBI 
show the fraction of our halo sample for which we found 
a main progenitor as a function of time. 

4.1. Density profiles during tidal mass loss 

The evolution of the mass distribution in satel- 
lite halos undergoing tidal stripping is often studied 
within an external f i xed potential (e.g. lDekel et al.ll2003l: 
lHavashi et al.ll2003t IKazantzidis et a l. 2004: Rea d et all 
120061 ) Theresolution tests in IKazantzi dis et all 1)20041 ) 
show that numerical effects lead to significant additional 
mass loss when an infalling subhalo is resolved with 
N=0.5xl0 6 particles and stripped down to a few thou- 
sand particles within a strong tidal field. The biggest 
subhalos in Via Lactea are almost as well r esolve d as the 
high resolution case in IKazantzidis et akl (|2004f ) . Many 
of the smaller ones lie far below their low resolution ex- 
ample and will suffer from artificial mass loss, especially 
when the tidal forces are strong, i.e. in the inner halo. 
Here we concentrate on the response to the tidal forces 
at pericenter passage of two of the largest, best resolved 
subhalos. 
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The first example is given is Figure 1101 This subhalo 
was accreted near a = 0.6 and completes three pericenter 
passages before z = 0. Its full track is shown in Figure 
[T2l while Figure [TOl depicts the mass distribution around 
this subhalo as it completes its second pericenter pas- 
sage at a = 0.844 (r por i = 7.0 kpc). The plot shows the 
mass enclosed within spherical windows of fixed physi- 
cal radii as a function of time. Not all of the enclosed 
mass will be bound to the subhalo. At the peak of the 
10 kpc line, for example, the majority of enclosed mass 
is associated to the underlying host, the density of which 
is is 5.5 x lOVrit = 8.3 x 10 6 M Q /kpc 3 at 7.0 kpc. On 
the other hand the host contribution to the mass en- 
closed within 1 kpc is negligible compared to the sub- 
halo's own contribution. The brief increase in enclosed 
mass for spheres with r ~ 1 kpc shortly after pericen- 
ter passage reflects a temporary contraction of the sub- 
halo as a response to th e rapidly varying potenti al, the 
so called tidal shock (e.g. lGnedm fe Ostrikerlll997f l. This 
contraction is only temporary, and shortly afterwards the 
mass in the affected spheres actually decreases. The en- 
ergy input from the tidal shock results in a net exp ansion 
(e.g. lHavashi et aT1l2003l : iFaltenbacher et al.ll2006h 9 . 

Particles in the outer regions of the satellite complete 
only a tiny fraction of their orbit around the subhalo 
center during the duration of the tidal shock. For these 
particles the tidal shock is impulsive and results in a 
maximal energy change. Particles near the subhalo cen- 
ter are less affected, since their internal orbital period 
is shorter than the dura tion of the shock. According to 
iGnedin fe Ostrikerll997l the energy input is proportional 
to 

AE(r) oc [1 + w(Y)t]~ 5/2 , (19) 

where t = 7rr P eri/V^eri is the duration of the tidal shock 
and ui = v c i TC (r) /2nr is the inverse of the circular orbit 
time in the subhalo. For the small pericenter in this ex- 
ample the shock duration is only r = tt (7.0 kpc) / (500 
km/s) = 42.9 Myr. This matches the subhalo orbital 
time at r oq = 0.20 kpc, i.e. at this scale cu(r)T equals 
one and AE(r) is reduced to 0.18 of the maximal value. 
Particles inside of r cq = 0.20 kpc should be less affected 
by the shock, unfortunately we cannot probe such small 
scales reliably. We do see however, that the shock is 
stronger in the outer halo: the mass within 1 kpc drops 
to a new constant value of 0.79 times the mass before the 
pericenter passage. Farther out the mass loss is larger: at 
a = 0.9 the mass within 10 kpc is only 52% of the value 
before the pericenter. Therefore the remnants of such 
a strong tidal interaction end up having lower densities 
at most radii, but with steeper, more concentrated den- 
sity profiles since more mass is removed from the outer 
regions. 

Before this pericenter passage the satellite has rymax = 
7.5 kpc. Its tidal radius at pericenter is much s maller, 
only 1.6 kpc. According to lHavashi et alj (|2003f ). satel- 
lites are fully disrupted when r t < 2r s ~ rvmax a t peri- 
center. This is clearly not the case in our example: the 
satellite survives this pericenter passage 10 , even though 

9 The expansion is large enough to overcome any tidal compres- 
sion suggested for subhalos orbiting in a $ oc r potential , as in the 
peer -1 inner region of the host halo (Dckcl et al. 2003). 

10 It even survives the subsequent, closer pericenter at only 5.2 
kpc. 
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Fig. 10. — Evolution of the mass distribution (solid lines) and 
tidal mass (dashed line) of a subhalo undergoing strong tidal forces. 
The thick solid lines show the mass enclosed within 1.0 and 10 kpc 
spheres around the subhalo center. The thin solid lines correspond 
to the mass within nine intermediate radii (1.3,1.6, 2.0, 2.5, 3.2, 
4.0, 5.0, 6.3 and 7.9). The smallest radius shown is 0.251 kpc, 
which is 2.8 times the force resolution. This halo approaches the 
Galactic center to within 8.3 kpc at a = 0.844 (vertical dotted 
line). Tidal mass loss is larger in the outer parts, but also the 
inner subhalo loses mass. In response to the strong tidal shock the 
satellite contracts just after the pericenter and expands soon after. 
The mass retained at a = 0.9 is 3.9 times larger than the mass 
within the tidal radius at the pericenter. A "delayed" tidal mass 
(dotted line) may be a better approximation to the bound mass. 
At r oq = 0.20 kpc (square) the subhalos internal orbital timcscale 
matches the duration of the tidal shock, (see main text for details). 

the tidal radius at pericenter is 4.7 times smaller than 
rymax- Figure [T2l shows that in general satellites survive 
even if they have several close pericenter passages like 
the one studied here. In ? 14.61 we find that total subhalo 
disruption happens very rarely, if at all. 

The second example is given is Figure [TT] This subhalo 
fell into the main host at a = 0.7 and completes only 
one relatively distant pericenter passage at a = 0.830 
(fperi = 58.3 kpc). Its track is shown in Figure [121 
The tidal forces at these distances are much weaker, and 
as a result there is no significant mass loss in the in- 
ner parts (less than 10% within 1 kpc) and only mild 
mass loss in the outer parts (29% within 10 kpc). A 
tidal shock can still be identified, but it is less strong 
and because of the longer shock duration it does not 
reach as far in as in the previous example. In this case 
t = 7r(56kpc)/(423km/s) = 406 Myr and r cq = 2.0 kpc. 
The increase at 10 kpc is purely due to the background 
density that peaks at 2.1 x 10 5 MQ/kpc 3 at pericenter. 
The shells around 3 kpc do show some contraction and 
expansion caused by the weak shock. Radii smaller than 
r cq = 2.0 kpc are practically unaffected by the shock. 
The mass retained at z = is larger (by a factor of 
1.8) than the mass within the tidal radius at pericen- 
ter: 10% of this mass is contributed by the host local 
background density at the position of the subhalo, and 
the remaining 90% can be attributed to the subhalo it- 
self. Remarkably, we find that the fraction of particles 
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Fig. 11. — Evolution of the mass distribution (solid lines) and 
tidal mass (dashed line) of a satellite undergoing weak tidal forces. 
This object falls into the host and has only one pericenter passage 
at 56 kpc and a = 0.8 3 (z = 0.20). The track of this subhalo is 
given in Figures 1121 and 1131 (dark green in the color version) . Tidal 
mass loss is large in the outer halo, while the inner subhalo remains 
unaffected. 

gravitationally bound to the subhalo (determined with 
SKID 11 ), is also 90%. Looking for additional bound ma- 
terial beyond the z = tidal radius increases the bound 
mass of the subhalo slighly: by 28% when going out to 
2.5 times the tidal radius. This larger bound mass is 16% 
higher than our tidal mass estimate, which includes the 
10% contribution from the host halo. 

The tests discussed above suggest that the simple tidal 
mass estimates are a good approximation (within 20%) 
to the bound mass at most times. During pericenter pas- 
sage our tidal mass may underestimate the bound mass 
by factors of about 2 to 4. The same problem would 
affect subhalo finders based on density only, since the 
bound material that is missed by the tidal criterion by 
definition lies near (within a factor of two) or below the 
host local background density. As a result subhalos at 
pericenter get assigned too little mass, and smaller sub- 
halos might be missed altogether. Our 6DFOF, in con- 
trast, always finds small groups, even when they lie below 
the background density. In the current implementation, 
however, we also do not assign enough mass to them. 

The transient dip in tidal mass during pericenter pas- 
sages seen in these two examples occurs because tidal 
stripping is not instantaneous; many particles remain 
bound, even though they lie beyond the tidal limit when 
the subhalo is near pericenter. Some semi-analytic sub- 
halo models incorporate this effect by removing only a 
certain fraction of the extra tidal subhalo mass at each 
time step St: 



Am d = M(>r t )6t/T s 



(20) 



T s is the time-scale for tidal stripping. We calculate the 
extra tidal mass M(> r t ) by subtracting mass that lies 



within the tidal radius in the current snapshot from the 
"delayed" tidal subhalo mass rrid at the previous snap- 
shot, rrid would be identical to the tidal mass in the 
limit of very rapid tidal stripping (T s — > St). T s is of- 
ten assumed to be equal to the satellite orbital time (e.g. 
iTavlor & BaLT[l[200ll : IZentner fc Bullockll2003[ ) 12 . Both 
of our (quite different) examples suggest that the strip- 
ping timescale is about six times shorter than the time it 
takes the satellite to complete one full orbit (T or bit). The 
delayed tidal masses assuming T s = T or bit/6 are shown 
with dotted lines in Figures [10] and [TT] This means that 
mass loss can be relatively quick, e.g. it is possible for a 
subhalo to lose more than half of its mass during only a 
tenth of the time it takes to complete one orbit. 

4.2. The hosts of Milky Way dwarf satellites 

A promising scenario to explain the low numbers of 
Local Group d warf galaxies relatiye to the abundanc e of 
CDM subhalos ([Moore et al.lll999l ; lKlvpin et al.lll999t) . is 
to suppress star formation in small halo s below a filter- 
ing m ass that increases after reionisation (IKravtsov et al 
20041) . Similar models (|Bullock et aLll2000l ; iMoore et al 



20061 ) select only systems above the atomic cooling mass 
at the reionisation epoch [z ~ 10). This too yields a re- 
alistic z = dwarf galaxy population and the disrupted 
building blocks are shown to match the spatial distri- 
butio n and kinematics of halo stars around t he Milky 
Way (jDiemand et all 120051 : IMoore et al.ll2006fl . In the 
Kravtsov et al. model the ten most luminous dwarfs are 
practically all found in the subhalos that had the largest 
peak circular velocities before these subhalos were af- 
fected by tides. The mild time-dependence of the filter- 
ing mass only leads to a few exceptions from this simple 
rule. 

To illustrate the evolution of the hosts of dwarf galax- 
ies in such scenarios we select two samples consisting 
of ten objects. The "largest before accretion" (LBA) 
sample is made up of the ten systems with the highest 
Knax throughout the entire simulation. These ten sys- 
tems all reached Vmax > 37.3 km/s at one time. The 
"early forming" (EF) sample consists of the ten halos 
with V,nax > 16.2 km/s (the atomic cool ing limit) at 
z = 9. 6. The EF sample corresponds to the IMoore et al.l 
(|2006| ) model, where sudden reionisation is assumed to 
have a strong effect on dwarf halos. The LBA scenario 
corresponds to allowing star formation only above a rel- 
atively high, constant critical size, a scenario of perma- 
nently inefficient galaxy formation in all smaller systems, 
independently of time-depe ndent changes in the e nviron- 
ment like reionisation. The IKravtsov et all (|2004l ) model 
would yield a selection that is intermediate between the 
LBA and EF samples. 

Six out of ten objects turn out to be the same in both 
samples. Because of the evolution of V mayi {z) during the 
growth of halos at high redshift, it is often the case that 
the ones reaching the largest sizes before accretion are 
also the ones that have the largest V max {z ~ 10). To 
avoid redundancy we plot the time evolution only for 
the EF sample, but show the z = properties of the 
LBA sample for comparison (Figures [T^] and [T3"|) . The 



11 Available at |http://v 
skid.html. 



-hpcc.astro.washington.edu/tools/ 



12 The time-scale used in IZentner et al.l (I2005T) is shorter: 
T'orbit/3.5. A factor of 2ir is missing in their Eq. 8 (A. Zentner, 
private communication) 
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Fig. 12. — Evolutionary tracks of the earliest forming substructure (EF sample, lines) and 2 = properties of the subhalos with the 
largest peak circular velocities before accretion (LB A sample, circles). Top l eft pa nel: distances to Galactic center versus cosmic expansion 
factor. Pericenter distances were calculated from the nearest snapshot (see £|4.5|| . Top right: evolution of tidal masses (or M200 f° r halos 
in low density environments, i.e. whenever M200 is smaller than the tidal mass). Bottom left: evolution of subhalo peak circular velocities. 
Bottom right: evolution of subhalo concentrations. Square with error bars gives the median and 68 % scatter of concentrations measured 
in field halos at z = 0. 

tween 1.5 and 4 times ^200 using all systems with V ma . x 
between 6.6 and 15 km/s). For systems with large mass 
loss, on the other hand, the concentrations increase with 
time and end up significantly above the field halo range. 
Figure [TJJ] illustrates the increase in concentration dur- 
ing tidal stripping in the V m ax — ry m ax plane. With our 
definition of concentration cy (eq. [7]), this parameter 
remains constant along V max oc rymax tracks. In the 
Vmax — fy m ax plane, halos start in the lower left corner 
at high rcdshift and then wander quickly towards the up- 
per right corner. During this active mass-growth phase 
7Vmax increases by a larger factor than Vmax, i.e. the con- 
centration cy decreases. After their active mass-growth 
phase they remain stationary in the upper right corner of 
this plane, until they experience tidal mass loss. Perhaps 
surprisingly, those satellites that do undergo tidal strip- 
ping retrace their paths in the Vmax — fvmax plane. This 
means that tidal stripping seems to exactly undo the in- 
side out subhalo assembly by removing mass from the 
outside in. Each stripped down z = remnant ends up 



subhalo histories are obviously very diverse. These ha- 
los have completed between zero and ten pericenter pas- 
sages. Some have lost over 99% of their mass, others 
no mass at all. The peak velocity V max may have been 
reduced by up to a factor of six, or not at all. Mass 
loss is often largest at the first pericenter passage (or 
the first few) and then becomes smaller except for the 
largest subhalo, the orbit of which is decaying quickl y 
because of dynamical friction (e.g. iTaffoni et all l2003h . 
All halos except one undergo pericenter passages that 
can l ead to tidally induc ed changes in galaxy morpholo- 
gies ([Maver et al.l l2001al fbl) . For dwarf host halos found 
near the Galactic center today tidal interactions were 
more violent, happened more often and started earlier, 
consistent with the very large mass-to-light ratios found 
for some inner dwarf satellites l|Maver et al.ll2007f ). 

Concentrations decline with time while halos are grow- 
ing. Later they remain constant for systems that lose no 
or only little mass and end up in the range where field 
halo concentrations are found (measured at z = be- 
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Fig. 13. — Evolution in the Vm ax — ?*Vmax plane of likely dwarf galaxy host halos (the EF sample, see text for details). Tracks start at 
high redshift when halos are still small, i.e. in the lower left corner of this plot. End points of the tracks at z = are marked with upward 
triangles. Only the end points are shown for the LB A sample (downward triangles). Each halo tends to grow and shrink within the same 
region of this plane, i.e. the end points resemble their (relatively concentrated) high redshift progenitors. Squares with error bars show the 
median location of z = field halos and the 68% scatter. Tidal stripping produces halos that are more concentrated than those in the field 
(i.e. smaller r Vmax at a given V m ax). 

resembling its own high redshift progenitor. Furthermore 
these z = remnants have high concentrations, typical 
of high redshift systems, and clearly higher than present- 
day halos that did not suffer significant mass loss. Both 
samples follow the concentration-radius relation (see Fig- 
ure [5]), i.e. the highest concentrations are found in sub- 
halos near the Galactic center. 

The differences between the two samples are small and 
presumably hard to detect from Local Group observa- 
tions: at z = the four EF halos that are not part of 
the LBA sample have smaller masses and lower Vmax but 
larger concentrations than the four LBA satellites that 
aren't part of the EF sample. Unfortunately, current Lo- 
cal Group dwarf galaxy mass models (e.g. lLokas et~aTl 
[20051 iKlevna et al] [20051 : iKoch et all [200% , based on 
radial velocity and surface bright ness data, can only 
place weak constraints on rvmax ([Strigari et al.l 120061 : 
iPenarrubia et all [2007V In both samples the hosts of 
the most luminous dwarf galaxies arc significantly more 
concentrated than field halos of similar mass: all ten 
cy's lie above the field halo median, and most (nine 
for EF, seven for LBA) lie above the mean halo scat- 
ter. Taking the higher concentrations of subhalos into 
account is importan t when estimating t he Vmax values of 
dwarf galaxy halos (IStrigari et al.ll2006l): using field halo 
concentrations instead ([Penarrubia et aLll2007l ) leads to 



higher V max values. 



4.3. Ensemble-averaged evolutionary tracks 

After the few individual examples discussed above we 
now go on to present ensemble-averaged evolutionary 
tracks, determined from a large number of subhalos. The 
tracks are selected by two criteria: the (sub)halo Vmax 
must be at least 5 km/s at some time, and the object 
must end up in one of the ten spherical shells around the 
main halo. 

Figure [TJ] show the median track of halos ending up 
within the first six shells, i.e. within r2oo = 389 kpc. Me- 
dian radii and 68% scatter illustrate where most subhalos 
lying in a given shell today were located in the past. Sub- 
halos in the two innermost shells at z — spend most 
of their time outside of these shells, i.e. most of these 
satellites are on more extended orbits and are near their 
pericenter at z = 0. The median radii in shells 3 to 6 are 
roughly constant since a = 0.5, i.e. the orbits of subhalos 
that lie within one of these shells today were distributed 
roughly symmetrically around this radial shell at earlier 
times. The scatter indicates that the typical peri- and 
apocenters lie outside the final shell, which is not sur- 
prising given the nearly isotropic orbits and the median 
ratio of apocentric to pericentric radii of 1:6 (see £ )4.5D . 

The radii also reveal some regularities between current 
positions and infall time. Halos in shell 1 today for ex- 
ample, are quite unlikely to have fallen in around a = 0.6 
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Fig. 14. — Evolution of (sub)halos that end up in the same radial shell at z = 0. Shells 1 to 6 are shown from top to bottom. Left panels 
show median and 68% scatter of (sub)halo distances from the center of the host halo. Due to their radial orbits, subhalos are often found 
outside of the fixed mass shell [shaded region) in which they end up at z=0. Right panels show median and 68% scatter of the fraction of 
current mass to final mass (dashed lines) and the one-fourth power of the fractions of peak circular velocity (solid lines). The thin dotted 
lines show the fraction of halos we are still able to trace at a given epoch. Halos without resolved main progenitor are not included in the 
median radii, but they are counted as zero when the medians of mass and peak circular velocity fractions are calculated. 



(and obviously also after a — 0.9). Those who did fall 
in around a = 0.6 are most likely found in shells 3 to 6 
at z = 0, after having completed one pericenter passage 
(for an example, see the track near 400 kpc at a = 0.6 in 
Figure 

In order to show how subhalo masses evolve, we fol- 
low two mass indicators over time. One is simply 
M(a)/M(a = 1), the ratio of the subhalo's mass at a 
given time to its mass today. The other is based on 
the ratio of peak circular velocities, [T4 nax (a)/K nax (a = 
l)] 4 . The later definition is motivated by the result 
that during tidal mass loss the peak circular velocity de- 
creases roughly like Kn a x cx M 1 / 4 (lHavashi et al.l 120031 ; 
iKazantzidis et all 12004 iKravtsov et al l l2004f ). When- 
ever this approximate scaling holds for our two mass in- 
dicators, they evolve proportional to each other. 

We have plotted the median and 68% scatter of these 
two mass indicators versus scale factor in the right hand 
panel of Figure [Ml In the outer regions near r2oo the two 
do indeed agree nicely, but closer to the halo center the 
masses are more strongly reduced than the peak circular 
velocities. At least in the innermost bin the tidal masses 
seem to underestimate the bound masses. From the me- 
dian radii (left panel) it is clear that most of the halos 
in this shell are at pericenter today, while in section [4.11 
we show that our definition of tidal mass tends to un- 



derestimate the true bound subhalos mass at pericenter. 
This explains the quick dip in the median mass fraction 
by almost a factor of 2 near 2 = 0. 

Not surprisingly, both the median mass and Vmax frac- 
tions show clearly that tidal stripping is stronger near the 
halo center. Stripping was also stronger at high redshifts: 
the median V max arc roughly constant after a = 0.7. In 
shells 3 to 6 both the median V max and median radii are 
roughly constant after a = 0.6. This supports the find- 
ings of Section [3j that subhalo and host halo evolution 
are closely linked. Early on the host system undergoes 
an active phase of merging and mass accretion, during 
which the subhalos arc accreted and their mass is re- 
duced quickly by tidal stripping at the first pericenter 
passage (Section I4.5j) . Once the host halo has formed 
there is little physical mass redistribution or accretion 
and the subhalo population becomes stationary. During 
this quiet epoch most subhalos move on stable orbits and 
their mass loss is relatively small. 

4.4. Formation histories and environment 

Figure [15] shows the ensemble-averaged orbital history 
of halos beyond r2oo today (in radial shells 7-10). Ac- 
cording to their current location these would be consid- 
ered "field" halos. However, many of them have actually 
orbited through the host halo at some earlier time. The 
resulting tidal interactions halted their growth and in 
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Fig. 15. — Same as Figure 1141 but for shells 7 to 10 that lie beyond r2oo = 398 kpc at z = 0. However, many of the "field" halos in shells 
7 and 8 were actually "subhalos" at some earlier time. The left panels contain include the fractions of halos that approached the galaxy 
halo to within 300 and 150 kpc anytime after z = 1. Such former subhalos have a reduced mass (and V ma x), due to tidal interactions 
with the main host at some earlier time. According to common definitions of halo formation time these stripped former subhalos would be 
classified as very early forming "field" halos. 

many cases even reduced their mass and 7 max - Such "for- 
mer subhalos" would be classified as very early forming 
field halos according to common definitions of formation 
time, because they reached a given fraction of their cur- 
rent, tidally reduced, mass or V max much earlier than the 
average "real" field halo that never passed through the 
host system. The population of "former subhalos" is sig- 
nificant: around galaxy clusters half of all halos found be- 
tween one and two virial radii today have passed through 
the c l uster at least on ce (|Balogh et al . 2000; Moore et al.l 
12004 iGill et al.ll2005[) . and for Via Lactea's subhalos this 
fraction is even larger, 0.74. The fraction is very similar 
when only the most massive (sub)halo orbits are used, 
i.e. we found no (sub)halo mass dependence. 

This illustrates clearly that a halo affects the forma- 
tion histories of many systems that are not within its 
virial radius (anymore). In other words, the assembly 
history of CDM halos must depend on their environment. 
Such correlations have indeed been quantified recently 
in ter ms of stronger cl ustering of early forming sub-M* 
halos (|Gao et al.ll2005| ) and as earlier median formation 
times in high density e nvironments (jSheth fc Tormenl 
12004 lHarker et al.ll2006h . Both studies are based on for- 
mation times defined relative to the z = halo proper- 
ties, analogous to our zg.5 definition (eq. [5]). This def- 
inition of formation time correlates strongly with mean 
environment density in and around the Via Lactea halo, 



as seen in Figure 1161 where we plot median formation 
times in our ten radial shells (cf. Figures [T4l and [T5|) as 
a function of the mean density in these shells. 

On the other hand, for a formation time based on pre- 
stripping halo properties, Zf orm (eq. [4]), this correlation 
disappears for halos outside today's virial radius. It is 
maintained inside the host halo because halos form be- 
fore they are accreted and median accretion redshifts are 
higher for subhalos that end up in the inner shells (sec 
Figure [T4]) . Note that the median Zf or m lies below z§5 
even in shells 9 and 10 where only very few halo (2.4 % 
and 0.006 %, respectively) interacted with the primary 
halo. This may be due to tidal interaction with a few 
other relatively large (Vmax — 70 km/s) halos in the out- 
skirts of Via Lactea. 

To summarize, we confirm th e (jBalogh et al.l 120001: 
iMoore et al.l[2004lGill et al.ll2005t) result that many sub- 
halos end up in the field (according to common defini- 
tions) and we illustrate that this leads to a clear envi- 
ronment dependence of halo assembly historie s like those 
found in (jGao et al.ll2"00l Barker et al.ll2006f ). A related 
explanation for the age dependence of halo clustering 
based on tid al interaction s with larger halos was recently 
proposed by IWang et alj l)2006f ). Defining halo forma- 
tion times relative to the maximal size a halo had over 
its lifetime removes the environment dependence of me- 
dian formation times. Nevertheless, this does not change 
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Fig. 16. — Median 2f orm (solid lines) and zgg (dashed lines), both 
with 68% scatter versus environment density. Using zg5, defined 
so that Vmax(z85) = 0.85V ma x(2 = 0), suggest that halos and sub- 
halo living in denser environments today (i.e. closer to the Galactic 
center) have formed earlier. The relation holds at all environments 
probed here, even well beyond the virial radius of the main halo. 
Using the pre-stripping halo size to define the formation time 2f orm 
the relation is only significant within the virial radius. 



the iGao et aTl (|2005f ) conclusion, that knowing only the 
z = mass of a halo is not sufficient to infer its accretion 
history or halo occupation distribution in a statistically 
correct way. Due to the strongly nonlinear tidal origin 
of this effect it seems difficult to correct the analytical 
approximations, suggesting that simulations should be 
employed whenever structure formation needs to be fol- 
lowed accurately. 

4.5. Mass loss per pericenter passage 

In this Section we quantify how much mass a subhalo 
loses per pericenter passage. We use the same sample of 
3883 relatively well resolved (Vmax > 5 km/s at some red- 
shift) subhalos as before. Pericenters are defined as local 
minima in the distance to the main progenitor. Only 
minima within iry max (z) are counted, to exclude min- 
ima caused by orbits around other progenitors. The time 
between stored snapshots (68.5 Myr) is too large to cap- 
ture all of the smaller pericenters. We calculate them by 
integrating orbits using the position and velocity of the 
subhalo and the spherically averaged mass distribution of 
the host halo at the nearest snapshot. The 68% interval 
of the directly measured r apo /r peY i distribution is nearly 
identical (within 0.01) as the corrected one, i.e. for most 
subhalos the calculated pericenter is similar to the one at 
the nearest snapshot. Our subhalos made up to 14 such 
pericenter passages, but most of them have completed 
only a few pericenters, or none at all (see Figure \T7\i . 
The ratios of pericenter and subsequent apocenter radii 
show that most subhalo orbits are quite radial. The me- 
dian ratio of our direct measureme nt (0.169) agrees ver y 
well with the derived value 1:6 from lGhigna et all (|1998f ). 
who used z = subhalo positions and velocities and 
the spherically averaged host density profile to integrate 
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subhalo orbits approximately The 90% interval extends 
from 0.035 to 0.666, i.e. only 5% of all subhalo orbits are 
rounder than about 2:3. These ratios are similar when we 
include only the last of several orbits (Table |4.5|) . When 
we include only the first pericenter passages the ratios 
are slightly smaller, the median is 0.133. The anisotropy 
parameter @(r) is about (3(r) ~ 0.55(r/r v i r ) 1 / 3 for these 
subhalos (and also for the total the dark matter) from 
r/r v i r ~ 0.2 to 1.0. The positive j3 values indicate radi- 
ally anisotropic subhalo (and dark matter) velocity dis- 
tributions, and the anisotropy increases with radius. 

Masses (and Kn ax ) are measured at the apocenter after 
the pericenter passage (t/), and at an earlier time so 
that the pericenter lies in the middle of these two times. 
This way we compare masses measured within similar, 
low background densities (even when a subhalo falls in 
for the first time) and we avoid the problem of underesti- 
mating the subhalo mass at pericenter (see N4.ip . We do 
not use late pericenter passages, i.e. when no subsequent 
apocenter is reached before z = 0. Mass loss is expressed 
as SM = [M(ti)-M(tf)]/M(U) and the reduction in sub- 
halo peak circular velocity 8V max is defined accordingly. 
Figure PT71 shows that the average mass loss (and the de- 
crease in Vmax) per pericenter is larger for orbits with 
smaller pericenter distances. The scatter is quite large, 
the rms scatter is about 0.22 in each radial bin (0.11 for 
6V max ). The average mass loss also depends strongly on 
the history of a subhalo: it is significantly larger when a 
subhalo passes through pericenter for the first time (Fig- 
ure [T7] and Table I4.5|) . The mass loss during the last of 
several orbits, on the other hand, lies significantly be- 
low the average over all orbits. Many of the individual 
tracks in Figure[T2l illustrate this behavior, i.e. they show 
a large early mass loss and nearly constant masses (and 
Vmax) near z = 0. The effect also manifests itself in the 
larger average mass loss before a = 0.6 in Figure [TH it is 
however smeared out because the first pericenters occur 
over a wide range of redshifts (the 68% interval extends 
from a = 0.31 to a = 0.63), which is earlier but overlap- 
ping largely with the distribution of all other pericenters 
(68% within a = 0.48 to a = 0.81). 

4.6. Tracing survival and merging forward in time 

The individual and ensemble-averaged tracks studied 
earlier by definition only include halos that have sur- 
vived (meaning they have a remnant above our resolu- 
tion limit) until today. This could potentially bias the 
reported median mass loss rates to lower levels. In this 
Section we quantify how many halos were stripped below 
our resolution limit and check if this reduces the mass loss 
reported for the survivors in the previous section. We se- 
lect halos with peak circular velocities above 10 km/s at 
z = 1 and look for their remnants today. Selecting only 
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Fig. 17. — Upper left panel: decrease in V max per orbit versus pericenter distance. Upper right panel: mass loss per orbit versus periccntcr 
distance. In both upper panels we show averages over all orbits (thick solid lines), over first orbits (upper thin lines) and over the last orbit 
of those subhalos that complete more than one orbit (lower thin lines). Pericenter distances are normalized to ?"Vmax( 2: peri)- Using a fixed 
scale like rymiif 2 = 0) = 69 kpc gives very similar results since rv ma x i s roughly constant since the last major merger at z ~ 1.7. The 
lower left panel gives the ratios of pericenter distance to distance at the subsequent apocenter and the lower right panel shows how many 
pericenter passages were completed by the 3883 subhalos in our sample. 



well resolved halos is necessary for two reasons: 

i) tidal stripping and destruction are overestimated 
due to numerical effects in bar ely resolved subha- 
los (the "overmerging p roblem" iMoore et al.lll996l ; 
iKazantzidis et al.ll2004f) . 



ii) one needs to identify the remnants even after severe 
tidal mass loss. 

At z = 1 there are 241 subhalos with Vi nax > 10 km/s 
within a sphere containing the final host halo mass (i.e. 
within shells 1-6). 232 of these are main progenitors of 
surviving subhalos, and two merge into a larger surviving 
subhalo between z = 1 and z = 0. Only the remaining 
seven subhalos are stripped below our resolution limit 
and disappear from our z — sample. Half of the de- 
bris from the tidally disrupted satellites are concentrated 
within a sphere of only 52 kpc around the Galactic cen- 
ter (for comparison, the half-mass radius of the halo is 
124 kpc), and the material beyond 52 kpc lies in three 
tidal streams oriented towards the center. Both the con- 
centration of the debris and the radial direction of the 
streams suggest that the destruction happened close to 
the Galactic center. Extending the sample size by in- 
cluding all 1542 subhalos with V max > 5 km/s at z = 1 
yields similar results: 2.4 % are lost and 1.3 % merge 
into a larger subhalo. 



Since about 97% of the subhalos selected at z = 1 sur- 
vive, the average evolutionary tracks of surviving systems 
given in Section 14.31 are representative for the majority 
of subhalos. It is also interesting to note that subhalo 
mergers are extremely rare, between z = 1 and z = the 
merger fraction is only about 1.3 %. 

Using the same subhalo selection we can also study 
what fraction of their z — 1 mass remains bound to their 
z = remnants. It turns out that this fraction depends 
strongly on the initial mass range of the selected subha- 
los. Larger subhalos retain less of their mass (Figure [TB")) . 
The most massive halo (light brown track in Figure IT!?)) 
has an orbit that decays quickly due to dynamical fric- 
tion and it loses most of is mass (98.9%) between z = 1 
and z = 0. Smaller subhalos are less affected by dynami- 
cal friction and lose significantly less mass. The increase 
in mass loss for subhalos with V ma , x (z = 1) < 7 km/s is 
likely artificial, and caused by insufficient numerical res- 
olution. Note that the z=0 subhalo velocity function of 
Via Lactea also starts to be affected by numerical limita- 
tions below the corresponding z=0 scale of about 5 km/s 
(Paper I). 

One consequence of larger mass loss in larger systems 
is that subhalo mass and velocity functions should be- 
come steeper with time, especially near when a region 
approaches virialisation, since this is the time when most 
of the tidal mass loss happens (Section [3]). Indeed, this 



Dicmand, Kuhlcn, & Madau 



19 



0.4 






0.35 


7 — — 


- 


0.3 


— - 


— - 


0.25 






s 

"~ 0.2 


- 


; 


0.15 






0.1 






0.05 










■- 



3 4 5 6 78 10 20 30 40 

V max (z=1)[km/s] 



Fig. 18. — Retained mass fraction from z = 1 to z = versus 
peak circular velocity range used to select the subhalos at z = 
1. The largest subhalo is represented by the square at 86 km/s. 
Subhalos were selected at z = 1 within a sphere containing Af200 = 
1.8 X 10 12 Mq. The largest subhalos are sinking towards the center 
due to dynamical friction, and they retain only a relatively small 
fraction of their initial mass. Smaller subhalos are on more stable 
orbits and lose less mass. The smallest systems (below about 7 
km/s or 1000 particles at z=l) suffer from an artificially large mass 
loss caused by the finite numerical resolution. 



trend can be observed in Figure [Jj The mass range used 
in this Figure extends down to Vmax = 5 km/s, the trend 
is stronger when only large subhalos arc considered: the 
slope of the cumulative velocity function of subhalos with 
Knax > 10 km/s measured within shells 1-6 grows from 
2.81 to 3.35 from z ~ 2 to z ~ 1, and remains roughly 
constant from then until z = 0. 

5. SUMMARY AND CONCLUSIONS 

We have analyzed the co-evolution of the Via Lactea 
host halo and its subhalo population. The simulation fol- 
lows the formation of a Milky Way-size halo in a WMAP 
3-year cosmology with 234 million DM particles. Here 
we summarize our main results: 

• In agreement with lPrada et alJ (|2006h we find that 
the formal virial radius, defined in terms of comov- 
ing scales, underestimates the actual virialized re- 
gion of ACDM galaxy halos. The increase in M200 
and -Mvi r after z = 1 is almost entirely due to ap- 
parent accretion, resulting from the artificial in- 
crease of the virial radius. Typically around 90% 
of the final M200 is already within the final r2oo at 
2 = 1. When halo mass is based on physical scales, 
such as Vmax or mass within rvmaxj we find no ev- 
idence for a late epoch of quiescent mass accretion 
as advocated by recent studies (e.g. IWechsler et al.l 
[20021: [Zhao et al.ll200l . 

• The collapse factors of shells enclosing a fixed mass 
are very different from the factor of two found 
in the idealized top-hat collapse. This causes the 
shortcomings of r v j r . 

• The abundance of substructure co-evolves closely 
with the host halo. The subhalo mass loss rate 



peaks between the epochs of turnaround and sta- 
bilization and declines after a region has virialized. 
Mass and velocity functions become slightly steeper 
during this process. 

• Tides remove subhalo mass from the outside in, 
which leads to higher concentrations compared to 
field halos of the same mass. This effect, combined 
with the earlier formation of inner subhalos, re- 
sults in strongly increasing subhalo concentrations 
towards the host center. 

• Selecting the earliest forming systems, or the 
largest before accretion, gives largely overlapping 
and at z = nearly indistinguishable subhalo sam- 
ples. They typically show large, early mass loss and 
high concentrations, especially those found near 
the Galactic center today. 

• We confirm the r es ult bv iBalogh et all (|2000t k 
iMoore et al.l ([2004D ; iGill et all ([20051 ) that many 
subhalos end up in the "field" (outside the virial 
radius) and quantify the environment dependence 
of halo formation times caused by this effect. Defin- 
ing halo formation times relative to the maximum 
circular velocity a halo reaches over its lifetime 
removes the environment dependence of median 
formation times, but not the environment depen- 
dence of halo mass assembly histories. Due to the 
strongly nonlinear tidal origin of the effect, cor- 
recting analytic approximations seems difficult and 
simulations should be employed whenever structure 
formation needs to be followed accurately. 

• At the first pericenter passage a larger average mass 
fraction is lost than during each one of the following 
orbits. The median peri- to apoc enter ratio is close 
to 1:6 (as in lGhigna et~allll998f ) and only 5 % of 
the subhalo orbits are rounder than 2:3. 

• We find that 97% of all z = 1 subhalos have a 
surviving z = remnant. The retained mass frac- 
tion is larger for subhalos with smaller initial mass. 
Satellites with Vmax — 10 km/s retain about 40% 
of their z = 1 mass at the present epoch. 
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